Non-Markovian stochastic description of quantum transport in photosynthetic systems 
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We analyze several aspects of the transport dynamics in the LH1-RC core of purple bacteria, which 
consists basically in a ring of antenna molecules that transport the energy into a target molecule, the 
reaction center, placed in the center of the ring. We show that the periodicity of the system plays an 
important role to explain the relevance of the initial state in the transport efficiency. This picture 
is modified, and the transport enhanced for any initial state, when considering that molecules have 
different energies, and when including their interaction with the environment. We study this last 
situation by using stochastic Schrodinger equations, both for Markovian and non-Markovian type 
of interactions. 
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INTRODUCTION 

Photosynthetic complexes are formed by ensembles of 
molecules that capture the energy from the sun, and 
through electron hopping transfer it to the so-called re- 
action center (RC): a particular part of the complex in 
which this energy is absorbed and used to produce chem- 
ical reactions [l| . 

Recent experiments have analyzed the so-called Fenna- 
Matthew-Olson (FMO) complex of the Green-Sulphur 
bacteria at cryogenic temperatures @, Q showing that, 
on the one hand, the transfer process between molecules 
occurs in a very short time scale (of the order of several 
hundreds of femtoseconds), and on the other hand, quan- 
tum coherence persist in the system during most part of 
the process. This resuls have been extended recently in 
[J], where it has been found evidence that quantum co- 
herence survives in FMO even at physiological tempera- 
ture for at least 300 fs. This has motivated a number of 
theoretical works that consider the photosynthetic com- 
plex as a quantum system, and try to analyse the basic 
mechanisms that explain the phenomena observed in the 
experiments. 

Most of these works analyse in particular the FMO 
complex [5|48|, and the LH1 — RC and LH2 complexes 
of the purple bacteria [9l4l2|, In both cases, molecules are 
subject to an environment of vibrating proteins, that can 
be described as a phonon bath. Since dissipation occurs 
in a much longer time scale, the bath produces basically 
dephasing [2|, [13J , a mechanism that gives rise to decoher- 
ence in the system and therefore leads to a recovering of 
the classical behavior. The question of how this dephas- 
ing affects the energy transport, and more generally the 
interplay between quantum coherence and dephasing in 
photosynthetic complexes has been widely discussed in 
the literature [Fj. Irl m, (14J16J . While as described above, 



some experimental results suggest that quantum coher- 
ences are present in the process, their role is yet to be 
determined. In some situations, like quantum walks in 
binary tree structures (see for instance Q and references 
therein) , quantum coherences might produce an exponen- 



tial speed-up with respect to the classical counterpart. 
However, when quantum coherences are too strong, they 
might give rise to a decrease of the quantum efficiency, 
i.e. the ability to transport the energy to the reaction 
center. Indeed, if the system is isolated from the envi- 
ronment and has uniform molecular energies, quantum 
destructive interferences may leave the system trapped 
in the so-called invariant subspaces |5|. Once there, it 
cannot scape and evolve into the target state in which 
all the energy is transferred to the reaction center. In 
this scenario, a certain quantity of environment noise, or 
alternatively some energy mismatches in the molecules, 
produces enough decoherence as to drive the system out 
of the invariant subspaces, leading to an increase of the 
transport efficiency [a, [g, |8[ . Environment assisted quan- 
tum transport in the FMO complex, as well as the entan- 
glement dynamics during the process have been recently 
analyzed for Markovian [7|, |l7| and for non-Markovian 
interactions [7|. 

In the present work we give a step further into these 
ideas, and by taking into account its symmetry proper- 
ties, we analyze the transport dynamics of the LH1-RC 
complex both for Markovian and non-Markovian cou- 
plings with the environment. Indeed, the LH1-RC core is 
a highly symmetric structure, as it was reveale d by pio- 
neering X-ray crystallography experiments (see [18[ for a 
review) . It consists in a ring of several antenna molecules 
coupled to a reaction center placed in the center of the 
ring Illll2|,ll9|,l20|. The periodic arrangement of antenna 
molecules allows us to describe their operators and their 
basis in the momentum space, a representation in which 
several interesting features can be observed. Particu- 
larly, we consider a simple case in which all the antenna 
molecules are coupled uniformly to a single reaction cen- 
ter. In this situation, it is shown that when the antenna 
molecules have uniform energies and are isolated from the 
environment, only the zero momentum component of the 
initial state is transferred into the reaction center. This 
component corresponds precisely to the symmetric part 
of the initial state, what explains qualitatively and con- 
firms some of the results obtained in [111, [l2| for the same 
system. In this scenario, the non-zero momentum com- 



ponents belong to an invariant subspace. However, when 
the molecular energies are non-uniform and/or the cou- 
pling with the environment is considered, it is shown that 
the energy contained in other momentum components 
can also be transferred to the reaction center, which leads 
to an improvement of the transport efficiency. We show 
also how some of these results can be extended to the 
case of several reaction center molecules non-uniformly 
coupled to the antenna molecules. 

The effects of the environment on the the transport 
dynamics is described by considering a quantum open 
system (QOS) approach. Within this framework, the dy- 
namics of the system is usually described with a master 



equation for its reduced density operator (see [2l| for a 
review), obtained after tracing out the environment de- 
grees of freedom. An alternative to this scheme are the 
so-called stochastic Schrodinger equations (SSEs). They 
evolve system vectors with the property that the average 
over their projector is equal to the density matrix. Some 
SSEs describe a deterministic evolution randomly inter- 
rupted by discontinuous quantum jumps, as described 
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in [22| for the Markov case and in 
Markovian case. The late equations have recently been 
applied to analyze the dynamics of the FMO complex 
[16J. A second type of stochastic equations lead to tra- 
jectories where the noise acts continuously at each time 
step of the evolution. This so-called quantum state dif- 
fusion scheme has been derived for both Markovian and 
non-Markovian interactions 
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, and it will be used 
here to describe some of the features of the LH1-RC com- 
plex. 

The stochastic scheme presents several advantages 
with respect to the use of the master equations: first, 
this formalism allows us to develop an approximation, 
the so-called post-Markov approximation, alternative to 
commonly used weak coupling approximation and valid 
for different parameter regimes. While the weak cou- 
pling approximation is valid for Tdeph <C 7, where 1/7 is 
the relaxation time of the environment, and l/T^eph is 
the system dephasing time, the post-Markov approxima- 
tion is valid for weakly coupled systems, such that the 
mean coupling J is such that J -C 7- Second, stochas- 
tic Schrodinger equations evolve the system state, that 
has say M components (M being the dimension of the 
system Hilbert space), while master equations evolve a 
matrix of M x M elements. Considering a single excita- 
tion in the system (an approximation that will be later 
explained in the context of photosynthetic systems), M 
is just the number of component particles or molecules. 
Hence, for a relatively large system like the LH1-RC com- 
plex, and provided that the number of stochastic trajec- 
tories needed to obtain the relevant quantities is not too 
large, SSEs might represent an advantage from the com- 
putational point of view. Note that, for instance, for 
M = 33 molecules and within the one excitation sector, 
a master equation deals with a Liouvillian of dimension 



33 4 ~ 10 6 , which is already very close to the computa- 
tional limit. 

The paper is organized as follows. First we analyze the 
system in the momentum space. Within this framework, 
the relevance of the initial states in the transfer efficiency 
is discussed. Then it is shown how the static disorder in 
the energies of the antenna molecules, as well as their 
coupling to the environment can increase this efficiency. 
In order to describe the system coupled to the environ- 
ment, SSEs are presented together with the post-Markov 
approximation. In that way the transport efficiency of 
the system undergoing dephasing is analyzed both for 
Markovian and non-Markovian interactions. 



A SIMPLE MODEL TO STUDY LH1-RC 

Let us now analyze the energy transport in the LH1- 
RC core of purple bacteria, based in a simplified picture 
of the more complex model described in [111, [12J and de- 
picted schematically in Figure (JTJ. For simplicity, the 
antenna molecules of the ring, as well as the reaction 
center, are considered as two-level systems. A more re- 
alistic situation in which the reaction center is modeled 

However, we expect 



by a molecular pair is found in [11 



that our simplified model captures some of the main as- 
pects of the quantum transport phenomena occurring in 
the system. On the other hand, the rate at which pur- 
ple molecules absorb photons (approx. one photon every 
two hours) is much smaller than the rate at which the 
energy is transfered to the reaction center. Thus, it is 
reasonable to assume that only one excitation is present 
in the system during the energy transport process. 

In this situation, the state of the system molecules 
can be described in terms of a basis {11)^} for j — 
1, • • • ,M + 1, where \l)j = |0 • • • lj • • ■ 0) is the state 
where only the jth molecule is in its excited electronic 
state . Note that antenna molecules are labeled with an 
index running from 1 to M, while the reaction center is 
labeled with the index M + 1. In that way, the effective 
system Hamiltonian can be written as 

M M 

h s = J2 uj i (7 t a 7 + J2 J ij a t a 7 +Yl r i( a t a M+i 

3=1 »>i=i j 

+ h.c.) + (wm+i + in)<TM +1 cr M+1 . (1) 

Here, u>i are the energies of the antenna molecules, while 
Jij = Jji are the electronic rates that describe the elec- 
tron hopping from site i to site j. Also, aj — |0)j(l| 
(similarly aj~) is the spin ladder operator describing the 
transition from one to zero (zero to one) excitations at 
site j, and cr^ +1 = |0)m+i(1| (similarly cr^ +1 ) is describ- 
ing the transition between one and zero (zero to one) ex- 
citations in the reaction center. The non-Hermitian term 
of the former Hamiltonian incf^ , i&m+i, has been added 



to phenomenologically describe the energy sink placed 
in the reaction center, which in addition has an energy 
ljm+i- The coupling coefficients of each molecule j with 
respect to the RC molecule are designed as Tj. 

As noted above, we consider here that there is a single 
molecule in the RC. However, our discussion can be easily 
extended to the case of having D > 1 molecules in the RC 
111 |l9| . In that situation, we should consider that the 
interaction Hamiltonian between each of the j antenna 
molecules and the each of the I molecules in the RC can 
be written as B L h\-rc = E;=i,d £j=i,m r '( cr / cr z7fic + 
h.c), where <Jj~ RC design the spin ladder operators corre- 
sponding to the RC molecules. This case will be analyzed 
in detail elsewhere. 

Let us now represent the Hamiltonian (jTJ) in the mo- 
mentum space. To this order, we consider the Fourier de- 
composition of spin operators for the antenna molecules, 
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(2) 



for j — 1, •••M. Here, rj = doj, and do the distance 
between different molecules in the ring. A similar expres- 
sion exists to relate the position and momentum basis, 

Having a system with periodic boundary conditions 
allows us to derive a crucial property, that may lead 
to important simplifications when expressing the system 
Hamiltonian in the momentum space. For a one dimen- 
sional periodic system like the molecules in the ring, we 
know that a J = crJ, M . According to ([2]), this means 
that e tqrM = 1, what implies that momentum q can only 
take discrete values q m = -j^jj, with m — 1, • • • , M. This 
leads to the equation 
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approximately aligned. In that situation, we can consider 
Ti sa T, so that T q w r<5 g o so that according to Q, only 
the zero momentum spin-wave, tro couples to the reac- 
tion center. This case corresponds to the configuration 
considered in the first part of the paper 11( , and is valid 



when the mean value of the set of Ti, is much smaller 



than their dispersion. 







M 



(3) 



FIG. 1: Schematic representation of the LH1-RC complex, 
with 32 antenna molecules distributed along a ring around 
the reaction center. The role of quantum coherence in the 
transport is intuitively explained as follows: When there is an 
excitation in a molecular ring with uniform energies, the quan- 
tum paths A and B interfere destructively with each other, 
so that path C (and therefore the transmission to the reac- 
tion center) is favored. Now, if we start with a symmetric 
state (corresponding to q = 0) this basic process will occur si- 
multaneously at every site, giving rise to a perfect transfer of 
the excitation to the reaction center. This destructive inter- 
ference is similar to the one that appears in quantum binary 
tree structures [8| . A schematic diagram of the basic building 
block of this process is shown at the right hand side of the 
figure. 



Using this property, together with the definition ([2J and 
a similar one for a^ , a simple form for Hs in the mo- 
mentum space can be obtained, 



(/.</' 



+ h.C.) + (WM+1 + iK ) (J M+l cr M+l- 



with u g>q > = J2< e 



i(i-i'). 
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3 = 1 



iqrj 



(4) 



(5) 



J,=EiJ: 



j"3< 



iqrj 



andT, = 4£,r 



q — m z-^j x o c 



iqrj 



Depending 
on the profile of T q , different momentum q are directly 
coupled to the RC. 

Let us assume an idealistic situation where all the 
dipole moments of the molecules within the system are 



On the one hand, if we assume uniform molecular ener- 
gies uij = ui, then U3 q>q i = uj8 q>q i and the Hamiltonian (0| 
does not describe any transfer process between different 
momentum components. In that case, if we start with 
an initial state that in the momentum basis is written as 
1^(0)) = £g- / 4g(0)|<z), the final transmission probability 
to the reaction center is Ptrans = -Pr(oo) = |„4 9= o(0)| 2 , 

withP(*) = i-ES 1 {4(*)° - 7(*)>- 

Hence, according to the picture described in 6], the 
subspacc spanned by the vectors {g}, with q ^ is an 
invariant subspace from which the excitation cannot es- 
cape and be transfered to the reaction center. This may 
explain in a very intuitive way some of the results of 
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12[ , in which it is shown that symmetrical initial con- 



ditions give rise to a better efficiency. Indeed, for a highly 
non-symmetric initial state, as |^(0)) = |1), that in the 



momentum basis has coefficients A q (0) 
transmission probability is just Ptrans = 



?*9 ri , the 



l/M. On the 



IM 



contrary, for a symmetric state as |\&o) 



IM 



E,K), 



with coefficients in the momentum basis A q = 5 q o, 
the transmission probability is the maximum possible 
Pr(t = oo) = 1. Figure [1] offers a possible explanation of 
the role of coherences in the process. 

On the other hand, if local site energies are 
non-uniform, the first term in the Hamiltonian Q, 
Eg «' UJ q,q' a q <J q~> leads to an energy transfer between dif- 
ferent momentum components. In that case, the q = 
component can eventually be repopulated by other mo- 
mentum components. This allows to inject more energy 
into the reaction center, and therefore to an increase in 
the transport efficiency. This can be observed in Fig. 
([2]) , where the evolution of the transmitted energy Pt is 
represented for uniform and non-uniform energies. 




FIG. 2: Evolution of the quantity P q =o(t) = (cr+ =0 (t)cr~ =0 (t)) 
and transmission probability, Pr(t), for two cases: A closed 
system in which every molecule has the same energy u (black 
and red solid lines for P q =o and Pt respectively), and a closed 
system in which molecular energies are randomly distributed 
as ojj = uio^j , with £j a randomly distributed number between 
and 1 and ojo = 20J (green and blue dashed lines for P q= o 
and Pt respectively). In the first case, P 9= o(0) gives the final 
transmission probability. In the second case the population 
of the zero momentum spin wave suffers some revivals pro- 
duced by the energy static disorder, in such a way that the 
transmission probability is increased. 

Throughout the paper, all the quantities are given in 
units of the typical transfer rate between near neighbours 
in the photosynthetic transport, J. According to |13| | typ- 
ical values for this quantity are J = 20cm~ 1 m 1.667ps. 

An alternative way to observe an increase of the trans- 
port efficiency is by realizing that molecules in the com- 
plex are not isolated, but rather they appear to be cou- 
pled to a phonemic bath produced by the surrounding 
proteins [29(. Since details about the microscopic model 
for this interaction are still unknown, we consider a phc- 
nomenological model Hamiltonian. As discussed in the 
introduction, the coupling produces mostly dephasing in 
the antenna molecules. Hence, our model Hamiltonian 



essentially accounts to modeling pure dephasing in the 
presence of a fluctuating field, and can be written as fol- 
lows 



M 





H 3 


= £ 

3=1 


A s Bj, 










(6) 
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and Bj = 


Ea h 


A (4 


+ 


ax), 


with 


a\ 



where Aj 

(at) the annihilation (creation) operators on the envi- 
ronment Hilbert space, and g^\ describe the coupling 
strength of the i molecule to the environment. The 
Hamiltonian © can be written in the momentum space 

as 



H } =Y. B ^ a \%' 



(7) 



where B q ^ = £\ ■e i{q ~ q ' ) B j . The Hamiltonian 
describes energy transfer between different momentum 
components of the antenna molecules, but only for inho- 
mogeneous couplings, when the environment operators 
Bj act differently over different molecules i. This means 
that the more uniform the interaction is, the less efficient 
is the mixing mechanism, what leads to a less efficient 
transport. 

As noted above, we have considered molecular dipolc 
moments are approximately aligned, which as analyzed 
in [ll|, [l9| , does not correspond to a realistic situation. 
However, we note that the representation in the momen- 
tum space here described, as well as some of the conclu- 
sions obtained, are useful even in the more general case in 
which the couplings V j (even for different reaction centers 
T l A are not uniformly distributed. To be more specific, 
we have introduced the data available in [19( for positions 
and dipole vectors of both antenna and reaction center 
molecules, in order to calculate the corresponding T l •. We 
have then computed each T l (as the Fourier transform of 
the set of T'- corresponding to each of the I — 4 reaction 
centers described in |l9|), and found that they are func- 
tions that present peaks at particular values of q. These 
peaks will correspond to the momentum values q that 
will be preferentially transfered to the reaction center. 
In addition, we have seen that for other values of q, all 
the functions T q are comparatively very small. Again, if 
there is no mechanism present in the system that mixes 
different momentum components (i.e. non-uniformities 
in the molecular energies, or dephasing), this momentum 
values in which r' sw will correspond to an invariant 
subspace: they will not be able transfer their energy to 
the sink. Further details of this realistic model will be 
analyzed elsewhere. 

In the next sections, let us analyse in detail the mixing 
mechanism produced by a dephasing, as well as its effects 
in the transport. 



REDUCED PROPAGATOR: A TOOL TO 
COMPUTE THE DYNAMIC OF A MANY-BODY 

QOS 

In order to describe the system dynamics, we assume a 
general model Hamiltonian H to t for a system with Hamil- 
tonian Hs, coupled to an environment with Hamiltonian 
Hb-, 



M 



Htot — Hs + Hb + Hj 

= H s + Y w > a A a > + Y Y 9j,x [ a l L ] + L ] a >fy) 

A j=l A 

where Lj is the coupling operator corresponding to the 
particle i and acts on the Hilbert space of the system, and 
the oj' x s are the frequencies of the harmonic oscillators. 
The Hamiltonian ^ is recaptured by considering Lj — 

The total wave function of the system, described with 
the Hamiltonian ©, evolves from its initial value | ^o) 
as, 



tf t )=W(t,0)| tf c 



(9) 



where U(t, 0) is the evolution operator in interaction pic- 
ture, 



U(t,0)=e ltiBt e 



iH B t —iH tot t 



(10) 



We can represent the state © in the basis of the envi- 
ronment, which is here chosen as the Bargmann coherent 
state basis [30(. In terms of this basis, ([9]) is written as 



**) = 



d(j,(zi)G(z± z \t) | -0o) I z\) 



(11) 



with the notation | zp) =| 2,9,1) | Z/3,2) ■ • • | z p,\) ' ' ' f° r 
the state of the environment, given by a tensor product 
of the states of all the A environmental oscillators. In 
our case, we are dealing with two states zq and Z\ corre- 
sponding respectively to /3 = 0, 1. In addition, we have 
defined the quantity dfj,(z) as a Gaussian measure given 
by 



Mz0) = ^VM 2 = 



n 



d 2 zp tX _ 



I^.aI 2 



(12) 



Equation (|llj) is obtained by introducing the closure rela- 
tion for the Bargmann coherent states, J dfi(z)\z)(z\ = I 
in ©, and considering as initial state \^f ) = \i J o)\zo), 
with \ip$) the initial state of the open system, and \zo) is 
the initial state of the environment. On the other hand, 
we have defined the function 26 28] 



G{z* 1 z Q \tO) = (z 1 \U{t,0)\z ), (13) 

that corresponds to the so-called system reduced propa- 
gator. This quantity acts on the system Hilbert space, 



giving the evolution of system state vectors from to t, 
conditioned that in the same time interval the environ- 
ment coordinates go from zq to z\. 

The phonon environment surrounding the antenna 
molecules in a photosynthctic complex is considered to 
be in a thermal state, and not in a pure state as \zq). 
Therefore, in order to describe the molecular dynamics a 
more general initial condition should be considered, 

Ptot(0) = / dn(zo)P T (zo,zZ)\zo){zo\® |Vo)(M (14) 

where Pt{zq, Zq) is the coherent state diagonal distribu- 
tion [30( corresponding to a thermal reservoir. 

Any expectation value of a system operator A can 
be obtained as (A(t)) = Trs(ps(t)A), with p s {t) = 
Tr B {pt t(t)) = Tr B {Jdfi(zo)PT(z 0l z*)U(t)\z )(zo\ ® 
\iJ)q){iPq\W (t). Introducing twice the closure relation of 
the environment, we get the following expression for p s (t) 
in terms of the reduced propagators, 



Ps(t) 



d(j,(zi) I / dfi,(zo)P(zo, z%) 

G(z 1 *z o |t0)|^ o )(^o|G- 1 (z^i|0i)) • (15) 



Note that this equation is just a representation of the 
reduced density matrix in terms of a coherent state basis 
of the environment. The reduced propagators are then 
quantities completely determined by the knowledge of the 
initial and the final state of the environment. However, 
this initial and final states are described by an infinite 
set of complex quantum numbers, corresponding to each 
of the harmonic oscillators that compose the environ- 
ment, distributed according to the measures dfi(zi) and 
dp J (zo)Pr respectively. Therefore, the integrals appearing 
in (|15[) are infinite and multidimensional, and a Monte- 
Carlo method should be used to solve them numerically. 
With this method, the integrals are sampled by choos- 
ing a random a set of coherent state coordinates z that 
are distributed according to their measure. In that situa- 
tion, the reduced propagator of the system G(z*,Zo\t,Q) 
become a stochastic object 26H28L |31|. 

We stress that the P function appearing in (fl5|) is by 
definition a quantum distribution, and hence it can take 
negative values. However, for a thermal environment it 
is positive definite, as it is required in order to perform 
the sampling over the coherent states. 

From equation (|15p . it is clear that two important in- 
gredients are needed in order to compute system quan- 
tum mean values with this method: first, we should be 
able to evolve the system propagator for different ini- 
tial and final states of the environment. In other words, 
for different noise histories. Second, we should perform 
the average over an ensemble of trajectories that is large 
enough to ensure that the Montecarlo integrals are sam- 
pled properly. With the Montecarlo method, the better 
the sampling is, the closer we are to the exact solution 



of the problem, except for the approximations made, in 
the equations of motion for the reduced propagators. We 
consider that the sampling is good enough when the re- 
sult converge to a value that remains fixed even if the 
number of trajectories included in the sampling is in- 
creased. 



L\z j0 , t ))G(zlz \tO)-T (20) 



where 



M r t 



T= V L\ f dTa pj (t-T)O j (r,t,zizo)G(z* 1 z \tO){21) 
„ „-_i jo 



p,j=i 



Evolution equation for the reduced propagator 

The evolution of the reduced propagator can be derived 
following a similar procedure as in 2JJ, |28j , and reads as 
follows (see further details of the derivation in Appendix 

A), 



dG{zlz \tO) _ 
dt 



L i z *j 



M 

= -iH s + J2 

3 = 1 

E , 
P .j=i Jo 

( Zl \u(t,o)L 3 (T,o)\z }, 



M - Llz j0t t) G(ztz \t0) 



L\ \ d,Ta P j (t - t) 



(16) 



where for each particle j, we have defined the functions 



z iP,t 

for (3 = 0, 1, and 

a pj (t - t) = 



A 



2^ 9p^9j,\t 
x 



-iuj\t 



-iui\(t-T) 



(17) 



(18) 



The quantity a P j (t — s) is the time autocorrelation func- 
tion of the noise Zji,t, as it can be easily verified by com- 
puting the average A4[z p itZ* lT ] with respect to the mea- 
sure dfj,(zi). From equation (|16p it can also be seen why 



the correlation function a P j is responsible of the depen- 
dence of the evolution of the system over its past history. 
Particularly, it is the kernel of an integral that goes from 
the initial time to the actual time t. In that situation, 
the slower the correlation function decays (i.e. the larger 
is the correlation time t c ), the more contributions appear 
from past times. 

In order to have a closed evolution equation 
for the reduced propagator, the matrix element 
(zi\l4(t,0)Lj(T,0)\zo) appearing in the last term of the 



equation (p~S|) should be expressed in terms of G(zlzo\t0). 
In order to do so, some approximation need to be consid- 
ered, such that the matrix element can be written as a 
certain system operator 0{z§z\ i7 t,T) multiplied by the 
propagator, 

(z\U(t,0)L 3 ( T ,0)\z ) = {z\Lj{T,ti)U{t,G)\*o) 

w 0,^,1, z* 1 z )G(z* 1 z \tQ).(l9) 

Inserting (fTTH) in (|16p . we get the following closed evolu- 
tion equation for the general reduced propagator, 



dG(z$z \tO) 
dt 



( 



M 

iH s + y](LjZ%., 



j'=i 



■ J 3 x 'ji,t 



Once an approximate form of Oj(t, t, z^Zo) is known, 
the equation (|2H)) can be used to integrate the reduced 
propagator along with its initial conditions G(z 1 zo|it) = 
exp(z^o). 

In the former equations, the operator Oj is constructed 

with the post-Markov approximation, first introduced in 

24] for a single particle quantum open system. This 

approximation is presented in more detail in the next 

section. 



POST-MARKOV APPROXIMATION 

We consider the post-Markov approximation in order 
to calculate a particular expression for Oj(r, t, z*Zq) in 
(fTTH) . To this purpose, we write the last term of (fTo]) as 

T = V / dTa PJ (t-T)Ll( Zl \L 3 (T,t)U(t,0)\z o m) 
p,j=i Jo 

and perform an expansion of Lj(t, t) = Lj(r—t) in Taylor 
series (t — r), 



L j( T -t) = L 3 



d 2 L 3 {r-t) 



dL 3 {r-t) 



dt 



(r-t) 



dt 2 



T=t 

2 , /n//_ +\3\ 



(r-ty + 0((r-tn (23) 



This expression, when inserted in (|22j) . gives rise to a 
series expansion of T of the form T — Y] n T^, for n — 
0, 1,2- ••, with 

p,j=i Jo 



where 



O r l(r,t\z* 1 z )G(z*z \t0)^(r~t) n 



{zi 



dL n Ar-t) 



dt n 



U(t,0)\z o ), 



(25) 



In the notation of (|19p . we can write Oj(t, t, z^zq) = 
J2 n Oj'( T > *> z i z o)- Note that inside the time integrals 
appearing in (|24|) . the maximum value for the quantity 
(t — t) n is the decaying time of the correlation function 
a p j(t), or correlation time r c . Hence, the n-th term of 
the series is of order r™ . 



The terms n = 0, 1 of the series for our system, de- 
scribed by the Hamiltonian ([1]), can be written as 

To=J2 dTa p] {t-T)LlL 3 G{zlz \m), (26) 

P,3 = 1 J ° 



with Lj = cfj<Tj, and 



A/ ,t 



71 = iJ2 I dTTll pj (T)G(ztz \tO) (27) 

pj=i Jo 



where 



K pj (t) = a pj (T)LlY^ Jjio-fai - Jijrfo 



(28) 



Within the single excitation sector and considering 
L p = cr^cTp, the former equation can be further sim- 
plified, since Lpa^ai = <r^ai5 P j. From (J26)) and (127)) 
we can infer that O®(rt\zlzo) = Q Art) = Lj, and 

Oj(rt\z*z ) = 0}{rt) = J2i( J ji (7 j'' Jl ~ J iJ a t a j)- Hence , 
for n = 0, 1 the operator Oj(t, i|z^zo) does not depend on 
the noises, and can be written as Oj(t, t) — Y^ n 0"(r, t). 
Further details of the calculation are explained in Ap- 
pendix B. 

In order to approximate T with a series expansion, it 
is important to establish the conditions under which this 
series converge, so that we can just keep the lower orders. 
This issue discussed in Appendix C. 



MASTER AND STOCHASTIC EQUATION FOR 

^ = L] 



Stochastic equations can be used to derive the cor- 
responding master equation [2J, [28(. This is done by 
considering the derivative of p s according to the defini- 
tion (|15p . then inserting the evolution equation for the 
reduced propagator (|16|) . and finally performing analyt- 
ically the averages over the environment degrees of free- 
dom. In that way, and for the simple case of an hermitian 



-t 



coupling Lj — ^,j — ^lj 
oder in n can be written as 



Aj, the master equation up to first 



dp s 

dt 



M . t 



i[H s , Ps}+Y] J d T a T p * (t - t) [A, , Ps O) (t, t) 
p,j=i Jo 

Y, / dTa^(t-T)[0 3 (t,T) Ps ,AA, (29 

p,j=i Jo 



with a pj (t) = J2\9j\9px(cosh(L)xP)cos(u)xt) - 
isin(uj\t)), and Oj(r,t) = Z) n =o,i < ^?( r ' i ) derived 
in the former section. Here, the constant j3 = 1/(KbT), 
with Kb the Boltzmann constant, and T the temperature 
of the reservoir. 



With some calculations, it can be seen that 
the master equation (|29[) . valid for an hermi- 



tian coupling, can be obtained also as -BzLl 



dt 



f t {J dp(z)G(z*0\tO)\^o)(ip \G Jf (z*0\t0)), by considering 
the following stochastic evolution, 



8G(z*0\t0) 



M 



dt 



iH s + J2 A jZ * t )G(z*0\tO) - 730) 



with 



M ft 

T~ V A, I dTu T p j(t - T)Oj(t,T)G(z*0\t0). (31) 



Note that this stochastic evolution is much simpler 
than ([2"U]). since it depends on a single noise such that 
f d/J,(z)z*z r — a T (t — t). For that reason, this will be 
the stochastic equation to use in this work. 



Evolution equations in the system basis 



We now write the stochastic (|30[) and the master equa- 
tion (|29|) in the system basis corresponding to the Hamil- 
tonian ([T]). 

In this basis, the evolution of the wave vector 
\ipt(z)) = S,-Oj(i)|l)j 5 where the coefficients Ctj(t) = 
(l\jG(z*0\tO)\ipo) represent the probability amplitude 
that the i molecule is excited, conditioned to an initial 
and final state of the environment and z respectively. 
For simplicity, we have excluded in the notation the de- 
pendency of aj(t) over these environment states. 

Considering the evolution equation for the propaga- 
tor G(z*Zo\tO) given by (|20|) . with the last term approx- 
imated with a post-Markov expansion up to second first 
order, T = X) n =o i Tni it is found that the evolution 
equation of a single coefficient a,j(t) for j = 1, • • • M is 



da 3 (t) 
dt 



M 

'I ^ Jj P a P {t) + WjOj(*) ) - iTa M+1 
+ z* t aj(t) - O jjaj (t) 

M 



" * £ ( ° X n J iP a pV) " °}v J pMt) ) > (32) 

and for the reaction center labeled as M + 1, 

da M +i(t) 



dt 



Here, we have defined 



-iluiM+l -^ JOAf+l(t) 

M 

* r £ a jW- 
3=1 



0™ = / dT{t-T) n C&(t-T). 

JO 



(33) 



(34) 



In this paper, we will consider the same correlation func- 
tion as in [13], defined as 



P3 



= S pjgj (-+i2j)e 7 * 



(35) 



where local noise and high temperature conditions, 
(3h-y = 0.25 < 1, have been assumed. This correlation 
function corresponds to an environment with a Drudc- 
Lorentz spectral density 



JAu:) =2g,- ? 7j 2 . 



(36) 



From equation (1291) . the corresponding master equation 
can be written as 



M 
' " " *X/ ( Jpl ^ ~ 5p,M+l)plj - JljO- - Sj,M+l)ppl 



dt 



i=i 



i(oj p - Uj)p p j - k(5 v< m+i + 5j,M+i))Pp 



M 



iT(5 p< M+l J^ Pi] + (1 - $p,M+l)pM+lj) 



1=1 

M 



+ iT{5j,M+\ /]ppi + (1 - 5j,M+\)PpM+i) 

i=i 

+ (o° p *-o°; + o PJ ~o° pp )p PJ 

M 



i=i 



M 



i J2 (o p ; - o)t + o}j - o} p ) j lpPlj , 



(37) 



i=i 



where p p j — (p\ps\j), and in this case {p,j} run from 1 to 
M + 1. Note however, that the last terms correspond to 
dephasing, and therefore only appear for {p,j} =/= M + 1. 



CHECKING THE STOCHASTIC EQUATIONS: A 
DIMER 



Before analysing the full LH1-RC system, we study the 
evolution of two molecules according to the Hamiltonian 
([1]) . We consider the evolution of the population of the 
molecule one, P\{t) = (cr|(t)o"i(i)), by using the master 
equation (|37|) . as well as the SSE (|32|) for different num- 
ber of trajectories NM. As shown in Fig. ^ when the 
number of trajectories NM is large enough both results 
are coincident. This illustrates the convergence of the 
Montecarlo method, as well as the fact that both schemes 
(master equation and SSEs) are basically equivalent. In 
other words, SSEs provide a good alternative to master 
equations to describe the dynamics of the system. 




FIG. 3: Comparison between master and SSE. The figure 
displays the time evolution of the population in the molecule 
1 of a dimer, Pi = {a[(t)ai(t)). Solid black line corresponds 
to the solution for the master equation, while red, blue, green 
and black dashed lines correspond to NM = 1, 100, 500, 1000 
trajectories of the SSE respectively. For both molecules, the 
coupling constant g = 0.3, the energies are wi = ui2 — 0, and 
the decaying rate of the correlation function is 7 = 10. The 
couplings are chosen as Jn = 1.5, J22 = 1 and J12 = 1.8. 
There is no reaction center in the dimer, so that T = 0. 



DEPHASING IN THE LH1-RC COMPLEX 

As noted above, when including dephasing in the de- 
scription of the transport process, a new mixing mecha- 
nism between different momentum components appears, 
what gives rise to an improvement of the transport effi- 
ciency with respect to the closed system case. In order to 
study qualitatively the effect of dephasing, we will con- 
sider that all the molecules have the same energy. This 
approximation is only valid when the molecular energy 
spread is small as compared to its mean value. However, 
it will allow us to describe the effects of the dephasing in 
the system dynamics in an separated way with respect to 
the effect of energy spread analysed in the former section. 

As we have seen in former sections that when con- 
sidering a closed system with zero dephasing, and uni- 
form energies and couplings to the RC, the total pop- 
ulation of non-symmetric momentum components q ^ 
0, PNs{t) — J2q=£o( a q (t)a q (t)} is a constant quantity 
through the evolution. This can be observed in the upper 
Fig. (J4|). On the contrary, when dephasing occurs, this 
quantity decreases severely, as can be seen in the curve 
corresponding to NM = 500 trajectories in the upper 
Fig. (H]) . This in turn is reflected as an increase of the 
population transfer Pt in an almost specular way (curve 
NM = 500 in the lower Fig. (]3|)). Along this section we 
will chose a smoothly decaying hopping, J p j — , _ 1 .. , , 

with do = 1/5, that phenomenology describes the de- 
caying of the hopping with the molecular distance. The 
coupling from antenna molecules to the reaction center is 
chosen as V — 0.5. As before, all quantities are in units 
of J. 

A brief comment concerning the sampling of SSEs is 



here in order. In both Figs. ((4]), it can be seen that 
for a number of trajectories NM > 250 the curves con- 
verge into a fixed value. This shows that for this number 
of trajectories, the Montecarlo integrals have been cor- 
rectly sampled, obtaining the final values for the relevant 
quantities. 




tn/Ji 



FIG. 4: Upper and lower figures represent respectively the 
evolution of the sum of the non-symmetric momentum com- 
ponents Pns and the transmission probability Pt for different 
number of stochastic trajectories NM. The inset of the later 
shows a detail of Pt for the closed system (g — 0) at short 
times. In both figures, solid black, red, blue and green cor- 
respond to NM — 10, 50, 250, 500 respectively, considering a 
Markovian situation, with F — 100, and a coupling constant 
g = 0.4. The dashed black corresponds to g — 0. 

We now analyze the total energy absorbed in the re- 
action center at Pr{t) at i = 5 for different decaying 
rates of the correlation function, 7 and for different cou- 
plings g. The time t = 5 is chosen in such a way that 
the system has already reached its steady state. In Fig. 
([5]) we analyze two situations: a large environment decay 
rate 7 = 100, such that for every coupling g displayed 
in the figure, the system is in the Markov regime, with 
7 ~S> Tdeph (Ta ep h being the system dephasing rate) ; and 
a situation in which non-Markovian effects begin to ap- 
pear, with 7 = 10. In this particular system, the results 
shown in Fig. (0 suggest that when dephasing is pro- 
duced by a non-Markovian reservoir, the transport effi- 
ciency is smaller than in the Markov case. The reason 
is that non-Markovian interactions preserve the coher- 
ence of the system longer, so that the environment is 
less effective in the task of taking the system out of the 
invariant state subspace, here spanned by the set of \q) 
with q 7^ 0. In other words, a non-Markovian reservoir 
is less efficient on mixing different momentum compo- 
nents than a Markovian reservoir. Although this behav- 
ior might be strongly model dependent (in the sense that 
different choices of the system parameters might give rise 
to different results), a similar effect is observed in [7| for 



the FMO complex. 

We should emphasize here that the post-Markov ap- 
proximation used in this work assumes that the electron 
hopping time scale is much larger than the correlation 
time of the environment. This condition might not be 
fulfilled in realistic photosynthetic complexes, and there- 
fore a more complete analysis is needed to describe the 
full non-Markovian behavior of the system. Some recent 
efforts in this direction can be found in FJ, I32H34 1 . 
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FIG. 5: Total energy absorbed in the reaction center Pr(t) 
at t — 5 for a Markovian case, 7 = 100 (solid black line), 
and a non-Markovian case, 7 = 10 (dashed blue). Here, we 
only reach intermediate values of the coupling parameter g, 
where the Zeno regime described in [g, |8| is not observed 
yet. We consider that this intermediate regime is the most 
intersting one, since it is the one where the transport efficiency 
is improved. 



CONCLUSIONS 

By considering a momentum space representation, we 
have analyzed in a systematic way several mechanisms 
that may affect the transport dynamics within the LH1- 
RC complex. To this order, we have considered a simple 
situation in which all antenna molecules are uniformly 
coupled to the reaction center, showing that, if addi- 
tionally the molecular energies are uniform and there is 
no dephasing, the final transmission probability corre- 
sponds to the initial population in the momentum q = 0, 
Pr{t = 00) = (o"g_ (0)<7g = o(0)). In other words, only the 
zero momentum component of the initial state is trans- 
mitted to the reaction center. On the other hand, we 
have shown how energy non-uniformities and dephasing 
produce energy transfer (or " mixing" ) from other mo- 
mentum q 7^ into the reaction center or sink, increasing 
the transmission efficiency. 

In order to analyze the effects of the dephasing in the 
system transport properties, SSEs are used to describe 
the system evolution. This scheme provides some com- 
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putational advantages with respect to master equation 
techniques, inasmuch as it evolves vectors instead of ma- 
trices. In addition to that, it allows to explore the system 
within the so-called post-Markov approximation, that is 
valid for environment correlation times shorter than the 
system electron hopping. When comparing the transmit- 
ted energy at t — 5[1/J] both for Markovian and non- 
Mar kovian interactions, it is found that at least in the 
present system, there is less transmission in the last case. 
This results are in accordance to what is found in [7J . In- 
deed, while dephasing destroys undesired coherences in 
the system, such that it can evolve out of the invari- 
ant subspaces and transmit more energy to the reaction 
center, non-Markovian interactions provide a mechanism 
to rebuild the coherences. Hence, non-Markovian effects 
may diminish the positive effects of the environment that 
exists for intermediate couplings. 

We note that even when the coupling of antenna 
molecules to the reaction center (or reaction centers) is 
non-uniform, we can still use the momentum represen- 
tation to determine which momentum values are most 
efficiently transfered into the reaction center, and which 
of them are not directly connected to it, and will need 
some momentum mixing mechanisms (like energy non- 
uniformities or dephasing) to transfer their energy into 
the sink. 

We thank D. Alonso, J.I. Cirac, A. Ekert, S. Huelga 
and M.B. Plenio for support and encouragement, and 
D. Alonso, M.C. Banuls, A. Chin, A. Ekert, G. Giedke, 
S. Huelga, M. Roncaglia, C. Navarrete-Benlloch and V. 
Vedral for interesting and useful dicussions. 



APPENDIX A: DERIVING THE EVOLUTION 
EQUATION OF THE REDUCED PROPAGATOR 

Let us give more details about the derivation of the 
evolution equation for the reduced propagator ([T5|) . 
From the definition of the reduced propagator, its deriva- 
tive is given by 



dG(zlz \tO) 
dt 



Z\ 



dU(t, 0) 



dt 



-'o 



(38) 



The evolution operator satisfies the Schrodinger equation 
in the partial interaction picture 



dU{t, 0) 
dt 



- iHs -i}^gj\(L]axe~ 
+L j a\e^ t )]u(t,0). 



iuj\tj 



(39) 



When inserted in (|55|) this equation leads to 

dG(z* lZo \tO) _ ( atj _ a ^ t n _^ xt . 



ot 



- iH s - i E L j9j,\ e 



l.A 



x G{z{z \tQ) ~ iL)Y,9^~ Wxt 

x (z\U(t,Q)a x (t,Q)\z ), (40) 

where we have used the property («i|al = (z\\z\ x . 
In the last term, we have replaced the matrix element 
(z\a x U(t,0)\zo) by (z\U(t,0)ax(t,0)\z ), with a A (t,0) = 
U~ 1 (t,0)a\L((t,0). This term can be further simplified 
by integrating the Heisenberg equations of motion for 
o A (t,0), 



dt 



ci\ 



(t,0) = -i^2 9j , x e- i ^ t L j (t,0), (41) 



with Ho = Hs + Hb , that leads to 

ft 



n\ 



(*,0) = a A (0,0)-iyV-,A / drL J (r,0)e^ r , (42) 

Jo 



with 



Lj(t',t) = e iHBt e -iH{t-t') L . e W{t-t') e -iH B t' _ (43) 

Plugging equation P2")l in (14TJ1) . we get the equation (fTBj) 
in the paper. 



APPENDIX B: SERIES EXPANSION OF T 

Let us now derive the terms T n with n = 0, 1, 2 appear- 
ing in the evolution equation (I2U1) . From the Heisenberg 
equation for Lj (r — t) , we can write 

< !MLZ3. = le lH ^ t e lH ^ T - t \H totl L 3 ]e~ lH ^ T ^e-' lHB (AA) 
dr 

Considering the a system Hamiltonian Hs given by (fT]), 
we find that 



J2Jij<rt(r-t)aJ(T-t)). (45) 



This quantity, evaluated in r = t as required in (|23|) can 
be written as 



dLj(T-t) 



dr 



= -i ( E J i*t°T E J ^t°j^) 



where we consider the notation a ■ (0) = ct . In the same 
way, the second derivative can be evaluated and leads to 



d 2 L 3 (T-t) 



-*E J y [(wj-w,-) + (B?-Bj)]^ 

i ^ 

/ j Jpl-Apj + 2 J Jpj-A-pi I (47 ) 
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a p a j 



<Jj~cr p , and Bj 



where we have defined A P j 

ExQjxiaxe-^^ + ale^). 

We now proceed to calculate the expressions for T n 
with n ~ 0, 1, 2, as defined in According the last results 
and equation (j24p . the zero order l~ n is simply 

% = J2 / dTa p] {t-T)LlL G{zlzu\m). (48) 
p,o=i Jq 

where G(z*Zo\tO) — (z\U(t, 0)|zo). The first order can be 
obtained by just inserting (|46|) in (|24|) , so that 



(49) 



Ti = iJ2 / dTTlZ pj (T)G(z^zo\tO) 

„ _-_i Jo 



Again, we may consider that (zi\U(t, 0)L/(t, 0)\zq) = 
(zi\Li(t, t)U(t,0)\zo), and a series expansion of the oper- 
ator Li(t, t) in terms of (r— t). Note that the term (|55l is 
a component of 72 , a second order term in the expansion. 
Hence, we just need to consider Li(r,t) ~ Li + 0(n) in 
(|51[) . The higher orders in n would give rise contributions 
of orders n > 2, which here are discarded. 

With this considerations, and replacing f[55]) in (fSTj) we 
get the following expression for the second order term of 
the expansion, 

72 = ( H / dTT 2 K p i{T)l ^2(gix - g P ,x)vx + (wi - w P ) 



P i 



o, 3 = l J0 



V f dTT 2 (V pk i{T) + Q p u(T)))G(z$z \tO) 
P ik Jo J 



(56) 



where in our case 



K pj {t) = ( J pj a P p(T)afa p - J P ja pj (T)a^aj 



where 



(50) 



Vx = «iV*" A * + zoxe-™^ +iJ2 L ' 



ie 



— ioj\t 



(57) 



The second line has been obtained after considering L*- = 
&t<7j and a single excitation in the system. The second 
order has the following form 

T 2 = V] / dTT 2 'R p i(T)lCi p {z* 1 z ) + (ui - Up)G(z*z \tO) 
„; Jo I 



APPENDIX C: VALIDITY CONDITIONS OF THE 
POST-MARKOV APPROXIMATION 



where 



and 



dTT 2 (V pk i(T) + Q pk i{T))G{z{z \tQ), 



C lp (zlz o ) = (z 1 \(B t l -B t p )U(t,0)\z o ) 



V p ki{r) = J p iJm {a P k(T)<T^ a v - a pp {r)a p a k } ; 
Q p h{t) = JpiJkp {a p k(r)a-t^r ~ a pi( T ) a i~ a k} ( 5S ) 

All the terms in the former equation are linear functions 
of the reduced propagator G(z*zo\t0), except for those 
containing the environment operators Bj, whose matrix 
elements should be correctly evaluated. In particular, we 
need to compute C p i{z\zq), 

Cpiizlza) = ^( 5i , A -.g PiA )(z 1 ^e^ t G(z 1 z |iO) 



In this paper we have analyzed the system dynamics 
by considering the so-called post-Markov approximation. 
(51) Let us now derive the validity conditions for this approxi- 
mation, and compare this conditions to the ones required 
in the weak-coupling approximation, and in the so-called 
Forster theory. 
(52) Indeed, the post-Markov approximation is valid only 

when the series expansion of T = Y] n T n converges at 
small n. This requires in particular that, for any j, the 
function 



Mt) = s ri 



dmm n a^Jm) 



(58) 



with S ~ ( Jjj , uij — lo p ) \raax is a monotonically decreasing 
function of n. The requirement can be written as 



F n {t) >F n+1 {t) 



(59) 
is the max- 



-iui x t 



zi\U(t,0)ax(t,0)\zo)) . (54) 



Replacing equation (J42I) in the matrix element 
(z\U(t,tO)a x (t,0)\z ), we get 



C p i(zlzo) 



for any n and any time t < t max , where t 
imum time at which the system is evolved. 

Let us now assume the environmental correlation func- 
tion (|35[) considered in this paper. With this correlation 
the function (1581 can be written as 



/■MH' f 11 " 1 ^ 11 



.*7] 



(60) 



Y / (9i.,x-9 P x)[Kxe t ^ t + z xe-^ t )G(z* 1 z \tO) 
a ^ where T[z] = L dte~ t t z ~ 1 is the Euler gamma function, 

iVe-^ t (z 1 |W(t,0)L ; (r,0)|z o )V (5%J d !>, *] = ^ dte'H^ 1 is the incomplete gamma 

; / function. Clearly, the function (I5U|) is a monotonically 
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decreasing function only if the ratio -<1. In terms of 
the electronic hopping J , the condition can be written as 



1 1 
J 7 



(61) 



In other words, it is required that the relaxation time of 
the environment is much shorter than the electron hop- 
ping. A further requirement is to chose 7 small enough, 
so that for the maximum time t max , the increasing func- 
tion in n, r[l + n] — T[l + n, £7], is compensated by the 

decreasing function I — J . From this discussion, it is 
clear that in principle within the post-Markov approx- 
imation the coupling between system and environment 
can be chosen arbitrarily. Therefore, it works even when 
the dephasing time scale of the system is comparable (or 
smaller) than the hopping time. In addition, for systems 
with small enough hopping, it may be even valid for cases 
in which the dephasing time of the system is comparable 
to the environment correlation time. 

The post-Markov approximation is an alternative to 
the weak coupling approximation and the so-called 
Forster theory, since it is valid in different parameter 
regimes. In the weak coupling approximation, the cou- 
pling strength between the system and the environment 
is considered to be small with respect to the magnitude of 
system and environment Hamiltonians (j < 5, B, where 
S = ||.ff,s|| and B = ||iJs||). Furthermore, along with the 
weak coupling approximation it is often assumed that the 
environment correlation time r c ~ I/7 is much smaller 
than the evolution time scale of the system undergoing 
dephasing l/^deph, i.e. 



1 1 

p »- 

-L deph 7 



(62) 



where Tg ~ g in terms of our parameters, so that the 
condition can be written as 7 S> g. In this approximation 
the coupling J and the energy detunings are arbitrary. 

In the Forster theory [13[ , a strong electronic coupling 
is considered, such that that electronic hopping occurs in 
a time scale much smaller than the evolution time scale 
of the system undergoing dephasing, 



1 1 

V > 7' 

1 deph J 



(63) 



or J » A. The correlation time of the environment, 
although somehow related with T$ is in principle irrele- 
vant. 
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